Coarsening Dynamics of an Antiferromagnetic model on the Kagome Lattice: 

Breakdown of the Critical Dynamic Scaling 
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We find a breakdown of the critical dynamic scaling in the coarsening dynamics of an antifer- 
romagnetic XY model on the kagome lattice when the system is quenched from disordered states 
into the Kosterlitz-Thouless (KT) phases at low temperatures. There exist multiple growing length 
scales: the length scales of the average separation between fractional vortices are found to be not 
proportional to the length scales of the quasi-ordered domains. They are instead related through a 
nontrivial power-law relation. The length scale of the quasi-ordered domains (as determined from 
optimal collapse of the correlation functions for the order parameter exp[3i0(r)]) does not follow a 
simple power law growth but exhibits an anomalous growth with time-dependent effective growth 
exponent. The breakdown of the critical dynamic scaling is accompanied by unusual relaxation 
dynamics in the decay of fractional (39) vortices, where the decay of the vortex numbers is charac- 
terized by an exponential function of logarithmic powers in time. 
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I. INTRODUCTION 

Thermodynamic systems quenched from a high- 
temperature disordered phase into a low-temperature or- 
dered phase exhibit characteristic growth of the length 
scale i of the ordered domains which in typical situa- 
tions can be represented as a power law £ ~ i^/^, where 
the growth exponent 1/z depends on the dimension of 
the space and of the relevant order parameter as well 
as on the conserved or nonconserved nature of the order 
parameter 

Phase ordering dynamics (or coarsening process) is 
usually accompanied by the annihilation and decay of 
the characteristic topological defects such as point vor- 
tices or domain walls which are generated in the initial 
disordered states. One of the most important notions in 
understanding and analyzing these coarsening processes 
is the so-called dynamic scaling hypothesis^ for the equal- 
time spatial correlation function of the order parameter, 
which is closely related to the observed self-similarity of 
the coarsening systems at different time instants. 

When the low temperature phase of the system is char- 
acterized by a quasi-long-range order with power law de- 
cay of the spatial correlation function of the order pa- 
rameter at equilibrium, then the dynamic scaling is gen- 
eralized to the critical dynamic scaling. One of the well 
known examples is the ferromagnetic Xl^modeP'^ on the 
square lattice. 

However, it should be noted that this (critical) dy- 
namic scaling hypothesis for the phase ordering (or quasi- 
ordering) dynamics has not been proved on some general 
theoretical basis. Therefore, it is not clearly known what 
is the condition for the validity of the dynamic scaling 
which is usually assumed in analyses of experimental re- 
sults or numerical simulations on coarsening dynamics. 



In typical situations, coarsening dynamics is investi- 
gated on systems where the phase ordering is accom- 
panied by the breaking of discrete or continuous global 
symmetry. However, there exist model systems that ex- 
hibit infinite ground state degeneracies of a discrete na- 
ture in addition to the usual (global) continuous sym- 
metry. Prominent examples are geometrically frustrated 
spin systems such as antiferromagnetic Ising models on a 
triangular lattice^i^, antiferromagnetic XYoi Heisenberg 
models on kagome^ or pyrochlore lattices.^ These sys- 
tems exhibit interesting equilibrium and noncquilibrium 
behavior due to the geometric frustration effect, includ- 
ing spin-glass-like relaxation dynamics without quenched 
disorder. 

In this work we perform numerical simulations on the 
coarsening of an antiferromagnetic XY model on the 
kagome lattice {KAFXY model)^'^ '^°i^^'^^ which is one 
of the simplest geometrically frustrated models with in- 
finite ground state degeneracies. Experimentally, this 
model can be realized in superconducting Josephson- 
junction arrays or superconducting wire networks ^^'^^i^^ 
on a kagome lattice when a perpendicular magnetic field 
of half fiux quantum (per plaquette) is applied on the 
system. We can also find physical examples in the 
anisotropic limit of Heisenberg antiferromagnets on the 
kagome lattice. It is well known that the system ex- 
hibits an infinite ground state degeneracy with finite 
entropjJ^. The system also exhibits a finite tempera- 
ture KT transition'^ corresponding to the unbinding of 
so-called fractional 30- vortices. That is, at low tempera- 
ture below the KT transition, equilibrium of the system 
will be characterized by quasi-long-range order of the or- 
der parameter ip^ = e^'^^i^) , Analogous to the case of 
simple ferromagnetic XY model on a square lattice, we 
might expect that the coarsening dynamics of the system 
would exhibit a critical dynamic scaling for the equal- 
time spatial correlation of the order parameter. 

Our simulations, however, show that the critical dy- 



2 



namic scaling is not obeyed very well at least for the 
time duration of our numerical simulations ranging sev- 
eral decades of time scale. That is, it was not possible to 
achieve a good scaling collapse for the equal-time spatial 
correlation functions of the 36 order parameter. Another 
signature of the breakdown of the critical dynamic scal- 
ing is that the length scales corresponding to the aver- 
age separation between fractional vortices are not propor- 
tional to the length scales of the growing quasi-ordered 
domains, which are instead related through a nontriv- 
ial power-law relation. This means that the two length 
scales exhibit different growth behavior in time. In terms 
of the decay of the fractional vortices, we also found that 
the fractional 30- vortices residing in the small triangular 
plaquettes exhibit faster decay, while, in contrast, those 
vortices sitting on the larger hexagonal plaquettes exhibit 
much slower decay that can be fitted by an exponential 
of logarithmic powers in time^''. 

In addition to these features of multiple length scales 
in the coarsening dynamics, the time dependence of the 
length scale of the quasi-ordered domains does not ex- 
hibit a simple power law growth but rather exhibits an 
anomalous growth with time-dependent effective growth 
exponent. This appears to be closely related to the un- 
usual relaxation dynamics in the non-power-law decay of 
the fractional {39) vortices. It is not clear yet whether 
the scaling may be restored in the asymptotic limit. 

II. THE MODEL SYSTEM AND SIMULATION 
METHODS 

The kagome lattice consists of corner-sharing trian- 
gles (Fig. [1]). In an antiferromagnetic XF model on the 
kagome lattice, the Hamilitonian is defined as 

H = - J^cos{ei~0j). (1) 

where J < 0, the sum runs over all nearest neighbor 
pairs of sites, and 6i denotes the angle of the planar spin 
at site i with respect to some fixed direction in the two 
dimensional spin space, (i, j) indicates all pairs of nearest 
neighbor sites in the kagome lattice. 

It is easy to see that the ground states of this system 
have the property that for all pairs of nearest neighbors 
i and j, the angle difference satisfies \9i — 9j\ = 2n/3 
which means that the sum of the three spins vectors on 
any unit triangle vanishes. Thus the space of the ground 
states of KAFXY model are equivalent (up to a global 
rotation) to the ground states of the three state Potts 
model on the kagome lattice. In Figs.[T](a)-(b) are shown 
two examples of simple ground states with long-range or- 
der, the so-called q = state and q = \/3 x \/3 state, 
where q refers to the wave vector corresponding to the 
periodicity of the chirality configuration. It is easy to see, 
however, that in addition to these ground states with sim- 
ple spatial order there also exist infinitely many ground 
states with no spatial order. It is well known that the 



system has a ground state entropy of So ~ 0.126fc_B per 
sitei^. Now, if we consider the angle variable 39 and the 
corresponding complex order parameter 'i>3 = exp{i39), 
the degenerate ground states are all completely ordered 
in terms of this new order parameter. And it has been 
shown that the system undergoes a KT transition at a 
finite temperature, where the spatial correlation of the 
order parameter exhibits an algebraic decay below the 
transition temperature. 

As for the KT transition temperature of this system, 
an analytic approximation was given as^^iii 

= J ^ 0.0756J. (2) 
72 ^ ' 

Numerical simulations were performed by Rzchowski^ 
where he found two slightly different estimates on the 
transition temperature, i.e., one estimate of Tc ~ 0.076 J 
based on the Binder cumulants of the order parameter 
and another of Tc ~ 0.070J from the helicity modulus 
(or correspondingly the decay exponent ?7(Tc) = 0.25 of 
the spatial correlation of the order parameter). 

In the present work, the coarsening dynamics of the 
model system is performed via kinetic Monte Carlo meth- 
ods with standard Metropolis algorithm for quenches to 
various temperatures near and below the KT transition. 
System sizes ranging up to x iV = 256 x 256 were em- 
ployed with periodic boundary conditions. In a kagome 
lattice, the number N of the total spins is 
The system is quenched from completely disordered ini- 
tial states down to a given low temperature with the pro- 
cess of coarsening being monitored through equal-time 
spatial correlation functions and the decay of topological 
defects, etc. 

In addition to quenches from disordered state to low 
temperatures, we also performed the so-called nonequi- 
librium relaxation by suddenly bringing the system from 
one of the ground states to some target temperatures 
around or below the KT transition. This method was 
found to be convenient for measuring the values of the 
critical exponent 77 for the equilibrium spatial correlation 
functions. 

When we let the systems evolve from random initial 
configurations the following quantities can be measured: 

1. The number density n„(t) of topological defects 
which are 36'- vortices in the KAFXY model 

...«) = (3) 

where Ny (t) is the total number of 30- vortices (both 
positive and negative) at time t and N denotes the 
total number of sites (i.e., spins). (•••) denotes 
an average over random initial configurations. We 
also count the separate number density of vortices 
residing on the hexagonal plaquettes (n/j) and those 
on the triangular plaquettes (rit). 
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2. The equal-time spatial correlation function of the 
36' order parameter -03(r,t) = e'^'^^'''*). 



C{r,t) = (^3*(r,i)^3(0,t)) (4) 




Y,eM^iO^{t)-3e,+r{t))) (5) 



3. Nonequilibrium spin autocorrelation functions 

= I ~ 0,{t))^ (6) 

Mt) = i(^exp(3z0,(O)-30,(t))\. (7) 



III. SIMULATION RESULTS 

We have performed dynamic Monte Carlo simulations 
of KAFXY model on a kagome lattice of dimensions 
256 X 256. In order to obtain the values of the equilib- 
rium critical exponents rj for the spatial decay of the equi- 
librium spatial correlation function Cgg {r) of the 30 or- 
der parameter, we have employed the so-called nonequi- 
librium relaxation (NER) method, where the system is 
suddenly brought from ground states to finite tempera- 
tures below or near the KT transition. Simulations were 
performed up to 655360 Monte Carlo steps which was 
sufficient for the equilibrium to be attained. This was 
confirmed by the collapse of the spatial correlations at 
later time stages. In Fig. [D^a) the correlation functions 
Ceq{r) are displayed on a log-log scale. The values of 
the exponent rj thus determined for temperatures rang- 
ing from T = 0.01 to T = 0.074 are shown in Fig. d^b) 
as well as in Table HI 
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V{T) 


0.01 


0.047(3) 


0.065 


0.254(11) 


0.02 


0.082(4) 


0.066 


0.259(13) 


0.03 


0.116(5) 


0.067 


0.264(12) 


0.04 


0.153(8) 


0.068 


0.268(12) 


0.05 


0.198(7) 


0.069 


0.280(11) 


0.06 


0.235(9) 


0.070 


0.300(13) 


0.061 


0.237(10) 


0.071 


0.308(12) 


0.062 


0.240(11) 


0.072 


0.310(13) 


0.063 


0.245(10) 


0.073 


0.330(15) 


0.064 


0.250(12) 


0.074 


0.340(13) 



TABLE I: The equilibrium exponent 77 for different tempera- 
tures for the power-law decay of the spatial correlation of the 
39 order parameter. 

As for the generation of the initial low energy con- 
figurations, we first select one of the three phases of 0, 



±27r/3 for each site at random and then applied Monte 
Carlo procedure at very low temperature of T = 0.005 
for 10000 MC steps with the restriction that, in each 
of the MC procedure, a new trial phase is chosen only 
among the above three phases. In this way, the final 
states obtained were very close (in energy) to a ground 
state. These states were employed as the initial states of 
NER in order to bring the system to the equilibrium at 
some given target temperature. 

From the graph of r]{T) shown in Fig. [H^b), we find 
that the exponent increases almost linearly in tempera- 
ture up to around T ~ 0.065, above which it starts to 
increase rather sharply. The linear regime can be fit ap- 
proximately by r]{T) ~ 3.9T. We also find that the expo- 
nent 77 takes the value of 1/4 near T{t] = 1/4) ~ 0.064 J 
(Table m. This temperature is expected to correspond to 
the ii'T transition, which apparently is a little lower than 
the theoretical^°'^^ or the numerical' estimates of earlier 
works mentioned above. 

Now, the coarsening dynamics of the model system 
under quench to low temperature from disordered ini- 
tial states is investigated by monitoring the equal-time 
spatial correlation functions of the tp^ order parameter. 
One of the most important features of typical coarsen- 
ing dynamics toward quasi-ordered phase is the critical 
dynamic scaling 

C(r,t)=r-"(^)/(r/L(t)), (8) 

where f{x) is the scaling function and L(t) is the growing 
length scale. It should be noted that the above scaling 
ansatz is based on the existence of a single growing length 
scale L{t). 

Figure [3]Ja)-(b) shows the equal-time correlation func- 
tions C{r,t) of KAFXY model for different time instants 
(from t = 10 to 655360) at the temperature T = 0.060. 
For these data, we attempted a critical dynamic scal- 
ing. The procedure of rescaling is as follows. For 
a given temperature, we take the value of ri{T) (de- 
termined as in Table |T| and evaluate the combination 
C{r,t) = r'^^'^^C{r,t). And then, for a given time in- 
stant, the length scale L{t) is determined by the condition 
C{r — L{t),t) = Co where the constant Co is chosen as 
Co = 0.2 or Co = 0.3 in this work. By plotting the scaled 
correlation functions C{r,t) in terms of the rescaled dis- 
tance r/L(t), we can check whether the critical dynamic 
scaling holds or not from the quality of the collapse of 
the scaled correlation functions. 

By this procedure, we found, rather unexpectedly, that 
the critical dynamic scaling does not hold in the coarsen- 
ing dynamics of KAFXY model, at least for the (Monte 
Carlo) time duration of our simulations (up to 655360 
MC steps) . The result of the scaling attempt is shown in 
Fig. ^c) for the case of T = 0.06 where we can see that 
the critical dyanmic scaling is not obeyed. We also tried 
different values of 77 for the rescaled correlation with no 
success. It might be possible that the scaling is restored 
in the limit of infinitely long time. 
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Even though the critical dynamic scaHng is not faith- 
fully obeyed, we may still extract approximate length 
scale L{t) of spatial correlation in the manner described 
above. The length scale L{t) thus obtained exhibits 
a rather unusual time-dependent behavior (Fig. [3l^d)). 
That is, in the early-time stage up to around t ^ 10^, L{t) 
exhibits a slow growth, which appears approximately in- 
dependent of the temperature. However, in the inter- 
mediate and late time stage of t > 10'^, the growth of 
the length scale L{t) becomes strongly dependent on the 
temperature. In addition, for a given temperature, no 
simple power law is found which is valid for the whole 
late time regime. Instead the local logarithmic slope ex- 
hibits a steady increase in time in the late time regime. 

We can investigate the time-dependence of the growth 
of the length scale by defining the effective local growth 
exponent as (3{t) = d\n{L{t)/d\n(t). Since we took the 
growing length scale L{t) only for discrete time instants 
ti with fixed interval in logarithmic scale, we evaluate the 
discrete version of the above logaritmic slope as 



\n{L{U+i)lL{U)) 



(9) 



where ti = 2% i = 1,2,- •• and = y/titi+i. Fig- 
ures. [3Ue)-(f) show the effective growth exponent at time 
t (for the case of Co = 0.2 (Fig. M.e)) and Co = 0.3 
(Fig.El^f)) respectively). We can see that, in the late time 
regime, the growth becomes definitely faster as the tem- 
perature is higher. Also, for the temperatures T > 0.05, 
the effective local growth exponent appears to be mono- 
tonically increasing in the late time stage. At higher 
temperatures, we can recognize some indications from 
Figs. [2Ie)-(f) that the effective slopes tend to converge 
around 0.5 ~ 0.6 in the latest time regime. It may be 
possible to interpret this as an indication of an asymp- 
totic behavior of diffusive growth. 

In order to understand the physical mechanism under- 
lying the breakdown of the critical dynamic scaling as 
well as the peculiar time-dependence of the length scale, 
we investigated the time dependence of the total number 
of the relevant topological defects, i.e., the fractional 36'- 
vortices residing on the triangular plaquettes as well as 
on the hexagonal plaquettes of the kagome lattice. The 
results are shown in Figs. ??(a)-(d) where we can see 
that, for a given temperature, the triangular vortices are 
decaying much faster than the hexagonal vortices. We 
also find that in general (for both types of the vortices) 
the decay of the total vortex numbers do not exhibit a 
simple power law behavior valid for the whole time range. 
One prominent feature in terms of the temperature de- 
pendence of the decay of the fractional vortex density is 
that, for the vortices on the triangular plaquettes, the 
decay rate increases as the temperature is lowered, while 
on the other hand, those vortices on the hexagonal pla- 
quettes exhibit opposite dependence on the temperature 
with the decay rate decreasing sharply as the tempera- 
ture is lowered. 



This can be interpreted as implying that there is al- 
most no (free) energy barrier for the motion and decay of 
the triangular vortices but that some finite barrier exists 
for the hexagonal vortices. We also find an interesting 
feature of the hexagonal vortices (at lower temperatures) 
where, at initial stage, they increase a little and then 
start to decay. This can probably be understood as the 
influence of the decay of the neighboring triangular vor- 
tices with their excess energy turned over to neighboring 
hexagonal vortices, thus generating hexagonal vortices. 

Therefore, we can conclude that, for wide range of tem- 
peratures, the relaxation of the fractional vortices resid- 
ing on the hexagonal plaquettes determines the coars- 
ening process such as the correlation length scale cor- 
responding to the quasi-ordered domains. In addition, 
the vortex relaxations exhibit a considerable deviation 
from power law behavior with slowly increasing (in abso- 
lute magnitude) local logarithmic slope in the late time 
regime. Especially in the case of the hexagonal vortices, 
the decay of the vortex number density could be fit to a 
form with 



Nh = Noexp[-b{\n{t)y 



Nq exp 



Vln(io) 



(10) 
(11) 



with b = In(io)^". Here, we found that typically a takes 
values around 3.3 ~ 4.0 (figure not shown here^*). We 
found that it could also be fitted to an stretched expo- 
nential form as 



Nhit)^NoeM~it/to)% 



(12) 



with a ~ 0.226 ±0.03. 

It appears that the breaking of the critical dynamic 
scaling is related to the anomalous relaxation of the num- 
ber of fractional vortices on the triangular and the hexag- 
onal plaquettes. In order to check this possibility, we an- 
alyzed the behavior of the length scales determined by 
the vortex distributions as follows. We define the length 
scale = l/yi9it which corresponds to the average sep- 
aration between the vortices on the trianular plaquettes. 
Similarly, we can define the length scale = 1/V^ 
corresponding to the average separation between vortices 
on the hexagonal plaquettes. We also define the length 
scales $,y = = l/^/uh + nt which represents the 

length scale corresponding to the total number density 
of vortices. 

Figure [5] shows the three length scales at T = 0.06 
which clearly shows the faster growth of the length scale 
for the triangular vortices compared with that for the 
hexagonal vortices. In Fig. ??(a)-(d) comparison is made 
between the length scales L{t) derived from the spatial 
correlation functions and the length scales ^« derived 
from the total vortex densities. We see that (for both 
cases of Co = 0.2 and Co = 0.3) there exists no simple 
proportional relationship between the two length scales. 
Rather, the length scale L{t) from the spatial correlation 
is seen to grow faster than the vortex lenth scale ^i, . This 
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probably implies that the vortices are not distributed 
evenly (statistically speaking) in the system but that the 
vortices are somehow distributed in a non-random man- 
ner such that some degree of clustering occurs. Interest- 
ing result is that the two length scales satisfy some non- 
trivial relationships such that £,v ~ L{t)'^ with A ~ 0.65 
for the case of Co = 0.2 (Fig. ??(b)) and A ~ 0.68 for the 
case of Co = 0.3 (Fig. ??(d)). This means that the vor- 
tex configuration of the system do not exhibit simple self- 
similarity at different time instants. Rather the vortices 
probably tend to cluster more unevenly as time passes by, 
leading to the correlation length scale L{t) growing faster 
than the length scale derived from the vortex density. 

The snapshot of defect configurations for the case of 
T = 0.06 is shown in Fig. ??(a)-(d). We can easily rec- 
ognize faster decay of those vortices on the triangular 
plaquettes. From the snapshots alone, however, it is not 
easy to detect the tendency of relative clustering of the 
vortices in the late time regime. 

Now we turn to the autocorrelation functions of the 
variables exp{i6) and exp(3?0). For these two variables, 
we find that the simulation results on the autocorre- 
lations are very different. For the case of exp(i36'), 
we obtain a approximately power-law decay behavior 
^3(0 ^ ^ in the early time regime up to around 
t ~ 10'^. The value of the exponent A ranges from 
A ~ 0.56 (for T = 0.01) to A ~ 0.90 (for T = 0.07). 
The exponent tends to increase slightly at higher tem- 
peratures. At later time, however, the autocorrelations 
decay faster than power-law relaxation of the early time 
stage. Due to the larger statistical fluctuations in the 
relaxation of the autocorrelation A3{t), we could not at- 
tempt a comparison of A^lt) with the growing length 
scale L{t). 

In contrast, in case of the phase exp(i6), the spin au- 
tocorrelaton exhibits less statistical fluctuations with a 
non-power law behavior— that can be reasonably fitted 



by 

A{t) = A„cxp[-b{\n{t)r] (13) 

with the exponent 7 — 2.0 ± 0.3 and 6 = In(to)^^. One 
of the simulation results is shown in Fig. ??(a-b) for the 
case of T = 0.06 with N, which shows a suitable fit to 
the above functional form with b ~ 0.15 and 7 ~ 2.0. 
Note that, in the limit of 7 = 1, A{t) reduces to a power 
law of A{t) ^ Aot^'^. The fitted values of 7 implies that 
A{t) exhibits a considerable deviation from a power-law 
behavior. We also attempted to fit the simulation results 
to stretched exponential form. This resulted in rather 
small values for the stretching exponent of a ~ 0.14. It 
is interesting to note that similar behavior of relaxation 
in the autocorrelation function of the order parameter 
was reported in the coarsening dynamics of the so-called 
Hamiltonian XY model on a square lattice.— 

In summary, we investigated the coarsening dynam- 
ics of the antiferromagnetic XY model on a kagome lat- 
tice. We found that the critical dynamic scaling is vio- 
lated. Novel and unusual coarsening dynamics may be 
attributable to the existence of infinite ground state de- 
generacies leading to nontrivial decay behavior of the 
density of defects residing on the triangular and hexago- 
nal plaquettes with nontrivial growth of multiple length 
scales. The competition between the critical thermal fluc- 
tuations of the equilibrium and the infinite ground state 
degeneracy probably gives rise to the nontrivial features 
of the relaxation dynamics. It may be possible to collapse 
the equal-time spatial correlation functions through a 
multiscaling scheme which we haven't tried yet. It would 
be also interesting to investigate the aging dynamics of 
this model system;^ 
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FIG. 2: (a) Equilibrium spatial correlation of exp(3i6) at various temperatures and (b) the corresponding exponents ri{T) 
vs. T, obtained from nonequilibrium relxation method. In (b), the dotted line represent 3.9T which fits reasonably well the 
behavior of r]{T) vs. T especially at low and intermediate temperature regime (T < 0.065). 
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FIG. 3: (a) Equal-time spatial correlation function C{r,t) for exp(3i^) for different time instants at T = 0.06, (b) the rescaled 
functions C{r,t) = r^{T)C{r,t), (c) a scaling attempt based on (b) with Co — 0.2, (d) growth of the length scale (using 
Co = 0.2), (e) the effective groYth exponents vs. time (using Co = 0.2), and (f) the effective growth exponents vs. time (with 
Co = 0.3). 
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FIG. 5: The length scales ^ft, and ^„ derived from the vortex densities versus time t for hexagonal, triangular plaquettes and 
the sum of these, respectively at T = 0.06. 
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